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' Recently, Gomez-Ullate et al. [1] have studied a particular Y-particle quan- 

go ' 

^) ■ field. They have shown that the Hamiltonian operator preserves a finite di- 

mensional space of functions and as such is quasi exactly solvable (QES). 
, In this paper we show that other types of invariant function spaces exist, 

which are in close relation to the algebraic properties of the elliptic functions. 
Accordingly, series of new algebraic eigenfunctions can be constructed. 
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I. INTRODUCTION 



The first example of a non-trivial, integrable quantum many-body Hamiltonian was found 
by Calogero [2]. It describes a system of N particles in one dimension interacting pairwise 
by means of an inverse square potential. The similar model endowed with an inverse sine- 
square potential is also integrable as shown by Sutherland [3]. In fact these two potentials are 
particular cases of a two parameter-family of potentials defined by the Weierstass function 
[4,5]. A detailed analysis of these models generalising the Calogero-Sutherland (CS) quantum 
models was reported in [6]. Their classical counterparts are discussed e.g. in [7]. 

While integrable models (classical or quantum) can be studied because of their mathe- 
matical interest, it became apparent in recent years that the CS models can be applied to a 
large number of fields of physics. These range from condensed matter (quantum Hall liquids, 
quantum spin chains, ..) [8] to gauge theories [9], soliton theory [10] as well as recently to 
questions related to black holes and (Anti)-deSitter space [11,12]. In particular it was shown 
in [12] that the asymptotic dynamics of 2-dimensional gravity in Anti-deSitter and deSitter 
space respectively can be described by a generalised two-body CS model. 

The property of a model to be integrable (i.e. to have a complete set of commuting 
constants of motion) does not necessarily imply that the spectrum and the eigenfunctions of 
the corresponding Hamiltonian can be constructed explicitely. The models which have this 
property are called solvable. From the beginning the CS models were known to be solvable, 
while further properties of their spectrum were obtained only recently, see e.g. [13,14]. How- 
ever, the explicit form of the spectrum is still missing as far as the full Weierstrass-function 
potential is considered for generic values of N. 

A step forward in the construction of solvable N-body problems interacting via a Weier- 
strass function was achieved in [1]. The authors indeed showed that, when the Weierstrass 
potential is supplemented by a suitable external potential, a finite number of eigenvectors 
can be computed explicitely in terms of special functions. Stated differently, the model is 
quasi exactly solvable (QES) according to the definition of [15]. In fact, the kind of inter- 
action considered in [1] and in the present paper generalises a potential first introduced in 
[16]- ^ 

Following the ideas of [15], the QES property holds when the Hamiltonian operator 
posseses a finite-dimensional invariant vector space of functions. Such a vector space was 
indeed constructed in [1] for the Hamiltonian considered. The purpose of this paper is to 
demonstrate that this Hamiltonian possesses alternative invariant finite-dimensional vector 
spaces of functions. The way these new vector spaces are constructed is very reminiscent to 
the multiple algebraisations of the Lame equations (see e.g. [17]), which occur due to the 
properties of the Jacobi elliptic functions. 

The Hamiltonian is presented in Sect. II. In this section we also give the transformation 
putting the Hamiltonian in a Lie-algebraic form which reveals its QES property. The new 
invariant vector spaces are constructed in Sect. Ill and the Hamiltonian is studied for 
particular values of the parameters. The results are summarized in Sect. IV. 

II. AN ELLIPTIC CALOGERO-SUTHERLAND MODEL 

The quantum Hamiltonian proposed recently by Gomez-Ullate et al. [1] is given by : 
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N d 2 

H N (x) = - ^ a - T + Viv(x) , x = (xi,x 2 , ....,xtv) . (1) 
It describes N particles on a line interacting through the potential 

N N N 

V N (x) = c m V(x k + if3) + 4b{b V ( 2x k) + a ( a ~ !) X) ^ + + " ^ • 

fc=l fc=l i,*=i 

(2) 

Here = V{z\g 2 ^gz) denotes the Weierstrass function with invariants g 2 , g%. The con- 
stants a, b are real and positive, c m is real. The term proportional to c m can be interpreted 
as the potential of an external field. 

The Hamiltonian (1) was shown to admit an invariant, finite dimensional vector space 
of functions [1]. Restricting the operator to this vector space, the eigenvalue equation 
Hpjip = Eip is reduced to a matrix equation and, accordingly, a finite number of eigenvectors 
can be determined algebraically. Following the definition of [15] the operator Hn is called 
Quasi Exactly Solvable (QES). 

To reveal this property has to be transformed appropriately. The authors of Ref. [1] 
introduced the function (called "gauge factor" ) 

Mx) = Winxj + 0) - V(x k + if3)f Y[[V'(x k + i(3)] b (3) 

j<k k 

and the new variables 

z k = V(x k + iP) , k = l,---,N. (4) 

Then a Hamiltonian Hn -spectrally equivalent to H^- is constructed according to Hn(z) = 
/x _1 (z)i/Af(x)/i(z). If the coupling constant c m is chosen according to 

c m = [2m + 2a(N - 1) + 46] [2m + 1 + 2a(N - 1) + 2b] , m e N (5) 

H N (z) preserves the finite dimensional polynomial space [18] 

N 

M rn = span{r[W^ ■ ■ ■ rft ; ^ ^ m > ( 6 ) 

i=i 

with the k-th elementary symmetric function 

r k = X z h z i2 ---z ik , l<k<N . (7) 

U<i2<---<H 

A lengthy calculation leads to 

N a2 N a I N a 

*«<■> = - E^ - *» E - (« + 2> + **<■> • («> 

fe=l k k,l=l K 1 K k=l K 

k^l 
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where p k = p(zk) and p' k = p'(zk) with 

p(z) = Az 3 - g 2 z - g 3 , p'(z) = £ = Mz 2 - g 2 . (9) 

In the following we will use the roots, e^, i — 1,2, 3, of p(z) : 

p(z) = 4{z - ei) (z - e 2 ) (z - e 3 ) = 4z 3 - - 9z ■ (10) 

These numbers are equal to the values of the Weierstrass function at its half-periods. 
The potential V N in (8) is given by : 

N 

V N (z) = m(12b + 8a( y N- 1) + 4m + 2)ri , Tl = ^z k . (11) 

k=i 

The crucial observation is that the Hamiltonian Hn(z) can be written as a quadratic poly- 
nomial of the differential operators 

d d N d 

Vk= g^' M jk = Tj—, U k = r k {r- J2 Ti ]fr. ) ' j,k = l,2,..,N (12) 

k k i 

with r = m. These operators form a representation of the Lie algebra sl(N + 1) for generic 
value of the real parameter r, for r = m they preserve the vector space M. m and the 
representation is finite dimensional. 

Denoting by H^ the part of Hn which increases the degree of elements of A4 m , we find: 

1 N 
= -4n(JV-m) (A/" + m + - + 2a(iV-l) + 36) , jV = ^ A4fc • (13) 

fc=i 

Obviously, the factor (A/" — m) leads to the annihilation of all the monomials in M. rn which 
have overall degree m. Therefore, we find: H^A4 m C M., m . As a consequence, eigenvectors 
of Hn (and therefore also of Hn) can be constructed in A4 m . In the following we will refer 
to this property as to an "algebraisation" of H N - 



III. ADDITIONAL GAUGE FACTORS 

Inspired by the construction of the Lame polynomials (see e.g. [17]), we introduce one 
further transformation of the Hamiltonian Hn '■ 

H N ^H N = jx^HNfL (14) 

with the gauge factor jl of the form : 

N 

fi(z) = H(z k - ei r(z k - e 2 y%z k - e 3 r • (15) 

k=l 
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The choice v\ = v 2 = vz = obviously corresponds to [1]. After a calculation, we find 
that for each value of the form 



Ui = or Vi = - - b , i = 1,2,3 (16) 

the Hamiltonian H N can be expressed as a quadratic combination of the operators (12) 
with suitable values (depending on the values of z/j's) of the parameter r. We then found 
eight gauge factors (15) leading to algebraisations of the initial operator H^. Let us now 
investigate the relations between r and the different parameters involved in the equations. 
We find that the degree- increasing part, say of is given by: 

H { + ] = -4n (V - (m + bn f - ^n f ^J (V + m + 2a(N - 1) + (3 - n f )b + ^(1 + n f ^j . 

(17) 

Here, rif denotes the number of non-zero exponents z/j, i — 1, 2, 3 in (15), i.e. is either 0, 1, 2 
or 3. Note that for rif = 1 and rtj — 2 three different algebraisations are available. 
If we allow m to be a non-integer and require instead that m with 

m = m + brif — -rif (18) 

is an integer, we conclude that now 

HnM^ C . (19) 

In the special case 6 = 0, we can distinguish two different cases: 1) both m and rh are integers 
and 2) only m is an integer. For 1) we find a quadruple algebraisation of the Hamiltonian 
H N (one algebraisation for rif = and three for rif = 2): 

HnMrn Q M m for u f — , (20a) 
H N M m -i C At m _i for n/ = 2 . (20b) 



Similarly, for 2) we find 



H N M m i C A^ m _i for n/ = 1 , (21a) 



^v-M m _| C ^ m _a for 77/ = 3 . (21b) 



Now, 777 = m — | should be an integer. Again, this is a quadruple algebraisation of the 
Hamiltonian H N (one algebraisation for 77/ = 3 and three for 77/ = 1). 
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A. a = b = : Relation between the Hamiltonian and the Lame operators 

In order to understand the pattern of the algebraic solutions obtained for the model (1), 
(2), it is useful to study the limit a = b = 0. Using the relation 

V(x + ip) = e 3 + (e 2 -e 3 )sn 2 (V^^x,k) , k 2 = ^—^1 ; (22) 

e i — e 3 

it is easy to see that for a = b = the operator (1) takes the form 

N d 2 

H N (u) = (ei - e 3 ) + 2m(2m + l)£; 2 sn 2 (^, fc)} + 2m(2m + l)e 3 N (23) 

3=1 i 

with Uj = \fe\ — e^Xj. The operator inside the brackets {} of (23) constitutes N decoupled 
copies of the Lame operator L{u) : 

d 2 

L{u) = -— + 2m{2m + l)k 2 sn 2 {u,k) , < k < 1 , (24) 

which admits (4m + 1) algebraic eigenvalues if m is an integer or a half integer. 

If m is an integer (m+ 1) eigenvectors of L(w) are of the form p m (sn 2 ) and (3m) eigenvec- 
tors are of the form cn-p m _i(sn 2 ), sn-p m _i(sn 2 ), dn-j9 m _i(sn 2 ). sn, cn, dn are abbreviations 
for the Jacobi elliptic functions sn(u, k), cn(u, k), dn(u, k) and p n denotes a polynomial of 
degree n in its argument. If m is a half integer 3(m +1/2) eigenvectors of L(u) are of the 
form sn • cn-j9 m+1 / 2 (sn 2 ), sn • dn-p m+1 / 2 (sn 2 ), cn • dn-p m+1 / 2 (sn 2 ) and (m — 1/2) eigenvectors 
are of the form sn ■ cn • dn • p m „!/ 2 (sn 2 ). 

Therefore, a total number of (4m + 1)^ algebraic eigenvectors of the Hamiltonian (23) 
can be constructed. However, not all of them are completely symmetric under the permu- 
tations of the coordinates. Since the procedure of algebraisation is crucially related to the 
symmetrized variables (see (7)), only the completely symmetric solutions can be hoped 
to be recovered in the generic case for which a^O and/or 6^0. 

Studying the solutions of the operator (23) and the structure of the eigenfunctions of the 
Lame operator, it is not difficult to see that the number of completely symmetric solutions 
is given by : 

Cm+N + (25a) 

if m is an integer and 

3C£, +JV + CZ +N -, ,m> = m+ l - (25b) 

if m is a half integer, respectively. C% denotes the usual combinatoric symbol. 

We find that for b = the number of algebraic solutions available by applying the method 
described here agrees nicely with these above numbers. Moreover, we checked for several 
particular cases that, indeed, the relevant Lame solutions are reproduced in the limit a — > 0. 
Note that in [1] only C^ +N solutions were found for integer values of m. Our supplementary 
factorisations therefore complete the pattern. 
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B. The case N = m = 2, 6 = 



For the choice N = m = 2, (20a) leads to a 6 x 6 matrix with respect to the basis 
{1,t 1 ,t 2 ,t?,t 1 t 2 ,t$} [1]: 

/ # 2 (2a + 26 + l) -2ag 3 Ag 3 \ 

16a + 246 + 20 £2(6+5) 4# 2 (a + 6+l) 2# 3 (1 - a) 

8a + 246 + 12 # 2 (2a + 26 + 5) -4# 3 (a + 1) 

8a +126 + 14 £2(6+2) 2g 3 

8a + 126 + 14 16(a + 36 + 3) # 2 (26 + 3) 

V 8a + 246 + 28 

(26) 

For (20b) we obtain 3 different 3x3 matrices with respect to the basis {1, Ti, t 2 } : 

/ (6 + 4a)e; g 2 (2a + 1) + 8ef -2a^ 3 \ 
a,= 14 + 8a (10 + 4a)e; £ 2 /2 + 4ef ,i = 1,2,3. (27) 
\ 28 + 8a (U + 4a)eiJ 

We thus obtain fifteen algebraic solutions, i.e. an additional nine to the ones obtained in [1]. 
In FIGs. la and lb we show the energy eigenvalues as functions of e for 

ei = 2 , e 2 = -l+e , e 3 = -1 - e (28) 

and a = and a = 5.0, respectively. FIG. la corresponds to two decoupled Lame operators. 
The limit e = further corresponds to the completely integrable case of two decoupled 
oscillators (e 2 = e 3 , so k = and the potential vanishes in (22), (23)). The eigenvalues of 
this system are of the form 3(jf + jf) — 40 where ji,j 2 are integers. The set of algebraic 
eigenvalues obtained with our factorisation (15) represents just the completely symmetric 
case, i.e. ji + j 2 = 2n, n = 0,1,2,... in this limit. This can be checked in FIG. la. In 
FIG. lb the effect of an interaction potential on the energy eigenvalues is demonstrated for 
a = 5.0. 

The case for which two of the numbers ei,e 2 ,e 3 are equal is in itself special, since the 
fifteen eigenvalues can be expressed as linear functions of a and the system is highly degen- 
erated, irrespectively of a. E.g. for e 1 = 2, e 2 = e 3 = — 1 three eigenvalues of the 6x6 
matrix are not degenerate : 

-8(5 + 4a) , -4(7 + 2a) , 8(1 + 2a) , (29a) 

the other three eigenvalues of the 6x6 matrix coincide with those of the 3x3 matrix hi : 

-8(2 + a) , 4(5 + 4a) , 8(7 + 2a) , (29b) 

and finally the eigenvalues of the 3x3 matrices hi and h 3 coincide and read : 

-2(17 + 10a) , -2(5 - 2a) , 2(7 + 2a) . (29c) 

This is clearly shown in FIG.s la and lb, where at e = three of the dotted curves, 
which correspond to three of the eigenvalues of the 6x6 matrix, and the three dashed 



7 



curves, which correspond to the three eigenvalues of hi, cross both for a = and a = 
5.0, respectively. Similarly, the three solid lines and the three dotted-dashed lines, which 
correspond to the three eigenvalues of the matrices h 2 and h 3 , respectively, cross at e = 0. 
How these degeneracies disappear for a generic choice of e i: % — 1, 2, 3 is also shown in these 
figures 

Finally, in FIG. 2 we demonstrate the dependence of the eigenvalues on the parameter 
a for the special choice ei = 2, e 2 = —3/2, = —1/2. 

IV. SUMMARY 

The construction of integrable models of Calogero-Sutherland (CS) type has recently 
received a lot of attention in relation to new applications related to different domains of 
theoretical physics. The class of N-body integrable models remains however very tiny and 
several generalisations are worth considering. The construction of quasi exactly solvable 
Hamiltonians describing N degrees of freedom appears to be a possible extension of the 
notion of integrable systems. As seen in [1,16] the potential can be more general than those 
related to the root system of a Lie algebra (typically of the type A N for potentials depending 
on the differences of the particles' coordinates). 

In this paper, we reconsidered such a QES model proposed recently in [1] . It depends on 
four parameters: two coupling constant a, b and the two periods of the Weierstrass function 
V, parametrized by gi,g%. More popular models are recovered for special limits of these 
constants: an Inozemtsev model for b = 0, a system of N decoupled Lame equations if 
a = b = and a system of N decoupled oscillators if, in addition, e 2 = e 3 (or equivalently 
g\ = 27 g%). We have seen that the case 6 = possesses a particularly rich algebraic spectrum. 

By investigation of the spectrum available in these limits, it appears that the solutions 
constructed in [1] do not constitute the full set of completely symmetric algebraic eigen- 
functions of the initial Hamiltonian (1). Following closely the construction of the Lame 
polynomials we have found additional algebraisations of the operator H N . The set of al- 
gebraic eigenfunctions obtained in this way coincides exactly with the number of possible 
algebraic functions. We assume that an extension of the type of Hamiltonian considered 
here to 2 x 2 matrix valued operators [19] might be possible, but leave this construction as 
a future project [20]. 

Note Added 

After the paper was finished several papers appeared dealing with the same topic. These 
are e.g. K. Takemura : math. Q A/0205274 and O. Chalykh et al. : math. QA/021 2029. 
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e 

FIG. la. The energy eigenvalues of the 6x6 matrix (dotted) and of the 3x3 matrices hi, 
(i = 1 dashed, i = 2 solid, i = 3 dotted-dashed), which correspond to the choice N = m = 2, are 
shown for a = b = as a function of e, where e\ = 2, e2 = — 1 + e and = — 1 — e. 
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e 

FIG. lb. Same as Fig. la, but for a = 5.0. 
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FIG. 2. The energy eigenvalues of the 6x6 matrix (dotted) and of the 3x3 matrices hi, (i = 1 
dashed, i = 2 solid, i = 3 dotted-dashed), which correspond to the choice N = m = 2, are shown 
for b = and ei = 2, = —3/2, e3 = —1/3 as a function of a. 
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